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Abstract 

A search for QCD-instanton-induced events in deep inelastic ep scattering has 
been performed with the ZEUS detector at the HERA coUider, using data corre- 
sponding to an integrated luminosity of 38pb~^. A kinematic range defined by 
cuts on the photon virtuality, > 120 GeV^, and on the Bjorken scaling vari- 
able, X > 10~^, has been investigated. The QCD-instanton induced events were 
modelled by the Monte Carlo generator QCDINS. A background-independent, 
conservative 95% confidence level upper limit for the instanton cross section of 
26 pb is obtained, to be compared with the theoretically expected value of 8.9 pb. 
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1 Introduction 



In the Standard Model, both the strong and the electroweak interactions are described 
by non-Abehan gauge theories. In such theories, the ground state has a rich topological 
structure, associated with non-perturbative fluctuations of the gauge fields P, called 
instantons [2]. They can be interpreted as tunneling processes of the gauge fields between 
topologically distinct types of vacuum states. 

Although the existence of instantons is required by the Standard Model, they have not 
been observed. While in electroweak interactions instantons are predicted to play a role 
only at centre-of-mass energies ^ lOTeV QCD instanton effects are expected to 

become sizeable at much lower energies, where they are predicted to have short-distance 
implications jSUlElIZI- In particular, they can induce characteristic events in deep inelastic 
scattering (DIS). The HERA ep collider offers a unique opportunity to discover instantons. 
Their discovery would constitute a confirmation of an essential non-perturbative Standard 
Model prediction, connected with the QCD vacuum. Results of such a search have recently 
been reported by the HI collaboration [Hj. 

Since instanton-induced events are predicted to contribute less than 1% to the cross sec- 
tion of the neutral current DIS sample used in this search, it is essential to find variables 
that efficiently discriminate between instanton and standard DIS events. Statistical dis- 
crimination methods have been employed to obtain event samples with a larger fraction 
of instanton events in a search for a possible signal. 

2 Experimental setup 

The data sample used in this search corresponds to an integrated luminosity of 38.3 ± 
0.6 pb~^ collected with the ZEUS detector at the HERA collider. During the years 
1996 and 1997, positrons of energy Eg, = 27.5 GeV collided with protons of energy 
Ep = 820 GeV. A detailed description of the ZEUS detector can be found elsewhere jH]. 
The main detector components used in the search presented here are the central tracking 
detector (CTD) ^U], operating in a magnetic field of 1.43 T provided by a thin supercon- 
ducting solenoid, and the uranium-scintillator calorimeter (CAL) [TT] . 

Tracking information is provided by the CTD, in which the momenta of tracks in the 
polar-angle^ region 15° < 6 < 164° are reconstructed. The CTD consists of 72 cylindrical 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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drift chamber layers, organised in nine superlayers. The relative transverse-momentum 
resolution for full-length tracks is cr{pT)/pT = 0.0058pT © 0.0065 © 0.0014/pr, with pt in 
GeV. 

The CAL covers 99.7% of the total solid angle. It is divided into three parts with a 
corresponding division in 6, as viewed from the nominal interaction point: forward (FCAL, 
2.6° < 9 < 36.7°), barrel (BCAL, 36.7° < 9 < 129.1°), and rear (RCAL, 129.1° < 9 < 
176.2°). Each of the CAL parts is subdivided into towers, which in turn are segmented 
longitudinally into one electromagnetic (EMC) and one (RCAL) or two (FCAL, BCAL) 
hadronic (HAC) sections. The smallest subdivision of the CAL is called a cell. Under 
test-beam conditions, the CAL single-particle relative energy resolution was a{E)/E = 
0.18/V^ for electrons and a{E)/E = 0.35/v^ for hadrons, with E in GeV. 

The luminosity was measured using the Bethe-Heitler reaction ep — e7p. The resulting 
small-angle energetic photons were measured by the luminosity monitor ^2], a lead- 
scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 

3 Characteristics of instanton-induced events 

Ringwald and Schrempp jUElEllZlElE] have identified kinematic regions in DIS that 
allow a perturbative calculation of instanton-induced processes. These processes lead to 
a characteristic final state, which may allow instanton-induced events to be distinguished 
from standard DIS processes. Figure ^ shows a diagram of an instanton-induced event 
in an ep collision. The incoming lepton emits a photon, with four-vector g, which in 
turn transforms into a quark-antiquark pair. One of these quarks hadronises to form the 
current jet with four-vector q" . The other quark, with four-vector g', fuses with a gluon 
(four-vector g) from the proton in the presence of an instanton. The phenomenological 
characteristics of instanton-induced events can be summarised as follows |3] : in the hard 
subprocess exactly one qq pair of each of the n/ kinematically accessible quark flavours 
participates in the quark gluon fusion process, either as incoming or outgoing fermion line. 
This gives rise to a high-multiplicity final state. The produced particles are expected to 
be isotropically distributed in their centre-of-mass frame. In addition, the events have 
large transverse energy in the hadronic centre-of-mass frame. 

Instanton-induced events in DIS were simulated using the Monte Carlo (MC) generator 
QCDINS 2.0 pi], which simulates the hard subprocess in the presence of an instanton. 
For the description of parton showers and hadronisation, HERWIG 5.9 jTHI is used. The 
simulation of the hard subprocess is accomplished by applying instanton perturbation 
theory around the one-instanton solution. By comparing the results obtained from in- 
stanton perturbation theory to those from lattice simulations of QCD, the fiducial region 
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of the QCDINS MC program has been derived [7j in terms of the kinematic variables 
Q'^ = —q'"^ = —{q — q")"^, x' = Q''^ / {2g ■ q') and = — g^, describing the hard subprocess 
of instanton-induced events: 

Q' > Q'min = 30.8 ■ Agl^ = VllSGeV^ ; x' > x'^^ = 0.35, 

where is the QCD scale. 

For QCDINS, cuts on the generated x' and Q'^ variables were made at the values given 
above in order to restrict the sample to the region where the calculation is reliable. 
Non-planar contributions [H], which are not taken into account in the calculation, were 
suppressed by a cut on the photon virtuality, Q"^ > Qmin — Q'min- 

The generated events were passed through the GEANT 3.13 -based ^7] ZEUS detector- 
and trigger-simulation programs. They were reconstructed and analysed using the same 
program chain as the data. 



4 Simulation of standard DIS events 



Standard DIS events were generated using the LEPTO 6.5 program ^H] interfaced to 
HERACLES ^ via DJANGOH 1.1 [201. The HERACLES program includes first-order 
electroweak radiative corrections. The CTEQ4 |25 proton parton distribution func- 
tions (PDF) were used. The QCD radiation was modelled with the colour-dipole model 
(CDM) [221 by using the ARIADNE 4.10 program |22]. Fragmentation into hadrons was 
performed using the Lund string model [21] as implemented in JETSET [22]. In order 
to improve the description of the sphericity distribution in the hadronic final state (see 
Section in)), the width of the transverse momentum distribution of primary hadrons (i.e. 
JETSET parameter PARJ(21)) was lowered to 0.28 GeV. The diffractive contribution to 
the neutral current sample was taken into account by adding 12% of diffractive RAPGAP 
events. This percentage was determined from a fit of the distribution of the variable 
^max |2E] , which is the pseudorapidity of the calorimeter energy deposit with the smallest 
polar angle and an energy above 400 MeV. In what follows, the term DJANGOH always 
refers to the combination of DJANGOH and RAPGAP. For a cross-check of the results, 
the generator HERWIG 5.9 [inj has been used; this program also provides the description 
of parton showers and hadronisation for the QCDINS MC. In HERWIG, coherence effects 
in the final-state cascade are included by angular ordering of successive parton emissions, 
and a clustering model is used for hadronisation j2Zj. Electroweak radiative effects are 
not included in HERWIG. Detector resolution and selection efficiency were simulated as 
for the signal sample (see Section El). 
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5 Event Selection and Reconstruction 



5.1 Reconstruction of kinematic variables 

Both track and calorimeter information were used for event reconstruction. Calorimeter 
cells were first grouped to form clusters which were then associated with tracks, where 
possible, to form energy-flow objects (EFOs) [2H]- The hadronic final state of an event 
comprises all EFOs that do not stem from the scattered positron. 

Positron identification was based on the pattern of energy deposits in the CAL [22] . The 
positron energy E^a was calculated using the double angle (DA) method [30] • Tracking 
information was used to determine the positron polar angle, if the polar angle in the CAL 
was above 0.3 rad and if the track traversed more than three CTD superlayers. Otherwise 
these angles were determined from the CAL information. 

The kinematic region investigated was defined by cuts on Q^, and on the Bjorken scaling 
variables x and y. The variables and x were reconstructed using the DA method, and 
y was reconstructed using the Jacquet-Blondel method p?T| . 

In order to reconstruct the kinematic variable Q'^ (see Fig. ^), EFOs were assigned to 
the current jet or to the instanton part of the hadronic final state. The current jet 
was identified by applying the /cT-cluster algorithm [32] in the longitudinally invariant 
mode jnSI on all EFOs in the hadronic centre-of-mass frame (hems). The photon direction 
was chosen as the negative Z direction. Monte Carlo studies showed, that in a region 
where the cross section for instanton related events is enhanced, and for high transverse jet 
momenta, the current quark has, on average, a smaller pseudorapidity than the partons 
assigned to the instanton. The current jet was therefore found as follows: a list of jets with 
a pseudorapidity t]^^^^ less than the transverse-momentum-weighted mean pseudorapidity 
of all EFOs in the hadronic final state was made {rj^^t^^ < XIefos V ' Pt/ XIefosPt)- Here 
the jet pseudorapidity was calculated from the jet four-momentum, obtained by summing 
the four-momenta of all the EFOs assigned to the jet. Of the jets in the list, the one with 
the highest transverse momentum was chosen as the current jet. 

The variable Q'^a "was calculated using the expression 



n'2 _ n2 I ^DA ~ Qba \^ p ^DA + Qda \^ „ ,^2 

QdA-QdA + ^^^ W^A 1^/^'^-^- 



jet 



with 



Wba = y Qda(1 - xba)/xba + , 
where Mjet is the current jet mass and nip the proton mass. According to MC studies, 
Q'^ is reconstructed with a relative accuracy of about 30%. However, the distribution of 
the reconstructed Q'^ has a long tail in the direction of overestimation of the true value. 
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5.2 Selection of neutral current DIS events 

A three-level trigger [HllSl] was used to select events online. A high-Q^ neutral current 
DIS sample was selected requiring, at the third level trigger, a positron with an energy 
greater than 4 GeV at a position outside a radius of 25 cm on the calorimeter. 

In order to select neutral current DIS events offline, the following cuts were applied, 
defining a fiducial sample: 

• kinematic cuts: Qda > 120 GeV^, xba > 10"^, yjB > 0.05; 

• cuts to ensure the quality of the positron reconstruction: 

o The radial position of the positron track impact point on the rear calorimeter 

surface was required to exceed 36 cm; 
o Eda > 10 GeV; 

o Ecnc < 5 GeV, where Ecnc is the energy, not associated with the positron, 
found inside a cone having an (77, 0) -radius R = y/(A0P~+7^^/P = 0.8 
around the track of the positron candidate; 

• suppression of photoproduction events: 

o < 0.90, where yd is calculated from the scattered electron; 

o 35 GeV < X]efos(-^ ~ Pz) < 65 GeV, where the sum runs over the energy- 

and Z-components of all EFO four-momenta; 
o DCA < 10 cm, where DCA is the distance of closest approach of the 
positron track to the centre of the cluster of GAL cells assigned to it; 

• vertex cut: Z position of the event vertex, |Zvtx| < 50 cm, consistent with an ep 

interaction; 

• restriction of the data sample to a region, where the QGDINS MG calculation is 

rehable: Qg^ > 140 GeVl 

Figure |21 shows the measured Q'^a distribution, compared to various MG predictions. 
There is agreement between data and the standard DJANGOH DIS MG sample at a level 
of about 10 %. A larger discrepancy is seen between the data and the HERWIG MG 
sample. Also shown is the prediction of the QGDINS MG program, normalised to the 
number of events in the data. It has a very different distribution from both the data and 
the standard DIS events. 

No cut was made on the variable x' for the standard DIS MG sample and the data because 
x' cannot be well reconstructed. On the other hand, the QGDINS Monte Garlo sample 
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includes a cut on the generated variable x' > 0.35. Lattice calculations show a steep 
decrease of the instanton contribution towards small values of x' < 0.35, corresponding 
to a small separation between instantons and anti-instantons, suggesting that this region 
can be neglected PH|I7|ITK]. 

Application of the above cuts resulted in a sample of 91846 events for the data. Normalised 
to the data luminosity, QCDINS predicts 578 events, DJANGOH predicts 88300 and 
HERWIG predicts 76400 events. The corresponding predicted instanton cross section 
is 8.9 pb. The statistical uncertainties on these numbers are negligible compared to the 
uncertainty on the luminosity. 

The discrepancy in the number of events predicted by HERWIG and DJANGOH can be 
traced to the cut in Q'^, and arises from differences in the Q'^ reconstruction due to the 
different hadronisation models in the MC programs. Without this cut, the numbers agree 
within the estimated uncertainty. 

In the subsequent analysis, numbers of standard DIS Monte Carlo events were normalised 
to the number of data events in the fiducial sample, where the predicted QCDINS contri- 
bution is negligible (~ 0.7%). 



6 Definition of discriminating variables 

Two kinds of discriminating variables have been considered: those connected with the 
kinematics of an instanton-induced event and those connected with the final-state particles 
of the instanton system, the so-called shape variables. The kinematic variables chosen 
were Q'da; (see Fig. ^ and Section E)) and pip", the transverse momentum of the current 
jet in the hems. 

The shape variables were calculated from a subset of hadronic final state EFOs assigned to 
the instanton, referred to as the "instanton region" in what follows. The instanton region 
comprises all EFOs which were not assigned to the current jet and lay in the hemisphere 
opposite to the outgoing proton remnant, i.e. polar angle 6^^^ > 90°. Once the instanton 
region was identified, the following shape variables were calculated: 

• A'^EFO, the multiplicity of EFOs; 

• A'^EFT; the multiplicity of tracks, that were used in constructing these EFOs; 

• C, the circularity; this is a measure of the isotropy of EFOs in the hems with respect to 
the photon-proton axis. To determine the circularity, the normalised two-dimensional 
momentum tensor, 

^a/ = y^7I^— ^ witha,/3 = X, F 
l^j \Pj,x + Pjx) 
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was computed in the hems from the three-momenta of the EFOs, requiring the num- 
ber of EFOs to be larger than two. From its eigenvalues Ai, A2, with Ai > A2, the 
circularity C was then obtained: C = 2(1 — Ai). 



S, the sphericity; this is a measure of how isotropically a collection of three-momenta is 
distributed. Large values correspond to a more isotropic distribution. The normalised 
momentum tensor, calculated from the EFOs in their centre-of-mass frame, 



.(3D) _ Yl,jPj,o'P'j,f^ 



MlT' = ^ , t-^^---^ ^ ^ with«,/3 = X,y,Z 



has eigenvalues Qi,Q2,Q3 with < Qi < Q2 < Q3. The sphericity is defined as S 
= 3/2 ■ {Qi + Q2)- The number of EFOs for this calculation had to be larger than two. 



• e', a measure of the density of the pseudorapidity 77 of EFOs in the hems; to calculate 
e', the EFOs were sorted with respect to the pseudorapidity r] in the hems, yielding 
A^EFO values ?7i, 772, • • • , Vnefo- The closed interval [r]i, rji^k], i + k < Nefo thus contains 
k + 1 EFOs. The variable e' is then defined by e' = ec — ^(A'efo — ^0) with b = 
0.339, No = 30 and 

k ^'^^-'^ 1 1 , f ^EFo/2 for A^EFO even, 

A^EFO - k ^ T]i+k -Vi I (A^EFO + l)/2 for iVEFo odd. 

A more detailed description and justification of the numbers b and Nq can be found 
elsewhere jSHl- 

Figures El and ID show the distributions of the discriminating variables in the fiducial 
sample (see Sectional). The contribution of instanton-induced events is 0.7 % according 
to the predictions of the QCDINS calculation. None of the variables have been corrected 
for detector or trigger effects. 

In general, there is a qualitative agreement between the shape of the data and the DIS 
MC samples DJANGOH and HERWIG. The QCDINS predictions, normalised to the 
number of data events, show quite different distributions, indicating the suitability of 
these variables for separating instanton induced events from background. 

However, the data are not reproduced in detail by the DIS MC simulations. Moreover, the 
two DIS MC descriptions differ from each other by a similar degree as from the data. For 
the DJANGOH sample, a parameter in the JETSET fragmentation program was tuned, 
as described in Section HI so that the sphericity distribution of the data is reproduced, as 
shown in Fig. |^. The description of the other variables, however, as shown in the figures, 
is not improved. 
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Given the uncertainties of the DIS MC predictions, and the smallness of the instanton 
contribution, it is not possible to make a rehable background subtraction. 

7 Instanton enhancement 

Several enhancement methods were compared to each other using the following criteria: 
From the numbers 

^o{Io) • number of standard DIS MC (QCDINS) events in the fiducial sample and 
Ne{Ie) '■ number of standard DIS MC (QCDINS) events in the instanton enhanced 
sample, 

the efficiencies r/, of QCDINS and r^r, of the standard DIS MC samples were obtained 
according to r/ = Ie/Iq and = Ne/Nq- Values of the separation power, Pg = rj/r^, 
were then compared for different samples of similar QCDINS efficiency. 

The enhancement methods investigated [HEI comprised a combination of one- dimensional 
cuts on the discriminating variables, an optimised choice of two-dimensional cuts and the 
Fisher algorithm [2Z1EH1; which performed best. 

The Fisher algorithm was used to separate standard DIS and instanton events by cutting 
on a linear combination of the variables Xj. For n input variables Xj, i = 1, 2, . . . ,n, 
the mean values a^** for the signal and xi'' for the background were determined. The 
correlation matrices between the variables were calculated according to 

^jfc.s = (a^i — Xi^){xl — Xk) 

events 

for the signal sample, and equivalently for the background {Vik^bj- These were averaged, 
Vik = |(^ifc,s + Vik,b), and the resulting matrix was inverted. The Fisher discriminant is 
defined as 

n 

t = ^WiXi , with Wi = )ik{xk' - Xk') ■ 

4=1 k 

Carli and Koblitz [21] have proposed a strategy to look for instanton-induced processes, 
using variables resulting from a careful optimisation procedure. An analysis with their 
variables, using the Fisher algorithm for signal enhancement, was also carried out. In 
addition, an analysis was carried out with a combination of their variables and those of 
Section IHl No significant improvement, as measured by the separation power P^, was 
achieved. 
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8 Background-independent limits 



In order to be independent of the standard DIS MC prediction, a conservative upper limit 
was set by assuming that all observed events after instanton-enhancement cuts are signal, 
i.e. the standard DIS background was zero. 

To derive limits on the cross section, the QCDINS MC program was assumed to give a 
correct description of instanton-induced events. For signal enhancement the Fisher al- 
gorithm was used, with the QCDINS sample as signal and the DJANGOH sample as 
background. In this procedure, the DJANGOH background sample is only used to deter- 
mine a good discriminating function; a non-optimal choice of background MC events will 
imply only that the instanton sensitivity of the measurement is not the highest achievable 
in principle. 

The Fisher discriminant was calculated from all six discriminating variables described 
in Section ini (5*, C, log{p^^^),NEFO, A^eft, e')- addition, a cut Q'da < 250 GeV^ was 
imposed, to improve instanton enhancement. 

The distribution of the Fisher discriminant t for the data, the QCDINS signal sample and 
the DJANGOH background sample is shown in Fig. |S| The distributions for the DIS MC 
simulation and the data are rather similar, with the curve of the instanton signal being 
well separated. The distribution normalised to the predicted instanton fraction is also 
shown. 

The numbers of data events above various values of the discriminating variable t are shown 
in Table ^ Also included are the numbers for the standard DIS MC samples. These are 
slightly larger than the number of data events. This is not surprising, since extremely 
restrictive cuts have been chosen, which only keep events in the tails of the standard DIS 
distributions. Problems are therefore expected in the modelling of the data. 

The number of data events kept by each cut, A^^^, was then compared with the theoretically 
predicted number of instanton events, Ath- The ratio A^/Ath is equal to the ratio of 
cross-sections R = ad/cth, where ath = 8.9 pb is the theoretically predicted instanton 
cross section for the cuts of the fiducial sample, and cr^ is the instanton cross section 
computed for the same cuts, assuming that all data events are signal. This is true under 
the assumption that, with zero background, the acceptances of the observed events and 
of the QCDINS MC sample are equal. Figure El shows the ratio R = crd/o"th, as a function 
of r/, the fraction of instanton-induced events remaining in the sample after various cuts 
in the Fisher t variable. 

In addition to the statistical uncertainties, uncertainties on the ratio R were taken into 
account by considering a ±3% change in the CAL energy scale, a change in the definition 
of the instanton region (see Section IH)), and using computed by the scattered electron 
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variables instead of Q^a- 

Another uncertainty stems from the cuts made on the variables Q'^ and x' by the QCDINS 
MC program at event generation, x' > 0.35 and Q'^ > 113 GeV^. These variables cannot 
be reconstructed well in the data. Therefore, the data might include some instanton- 
induced events with real x', Q'^ values below these cut values. If, with the help of an 
accurate reconstruction of Q'^ and x', these events could be removed from the data sample, 
it would lower the ratio R, yielding a value closer to the theoretical prediction. 

From the value of R and its uncertainties, an upper limit for the instanton cross section 
can be derived for any given choice of r/. For example, at r/ = 10%, the upper limit 
would be 30 pb, to be compared with the theoretically predicted value of 8.9 pb. In a 
similar analysis by the HI collaboration |H], in a region of phase space with a lower range 
of g2 (lOGeV^ < < 100 GeV^), and at a comparable value r/ = 10%, an upper limit 
of 221 pb at a 95 % confidence level (c.l.) was reported. This value is about a factor of 
five above the corresponding theoretical prediction. 

In order to derive a conservative upper limit without an explicit choice of r/, the t distri- 
bution was computed for a range of values of the instanton cross section. For a specific 
choice of the cross section, the t distributions for data and instanton events become equal 
at a certain specific value to, the instanton distribution overshooting the data for t > t^. 
An upper limit can then be set by choosing the instanton cross section such that the 
number of instanton events exceeds the number of data events for t > ai a, (one sided) 
95% c.l. This method yields an upper limit at a 95% c.l. for the instanton cross section, 

CTinst, of 

ffinst < 26 pb, 

to be compared with the theoretically predicted cross section of 8.9 pb. The fraction of 
instanton-induced events remaining after the cut was = 4.6%. 

9 Conclusion 

A search for QCD-instanton-induced events has been performed in neutral current deep 
inelastic ep scattering, based on an integrated luminosity of 38pb~^ in the kinematic 
range > 120 GeV^, x > 10~^. Cuts on the Fisher discriminant have been used to 
obtain instanton-enhanced subsamples. 

Assuming that all data events belong to an instanton signal, a conservative background- 
independent upper limit on the instanton cross section of 26 pb at a 95% c.l. has been 
set, to be compared to the theoretically predicted cross section of 8.9 pb. 
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ri [%] 


DATA 


QCDINS 


DJANGOH 


Ps 


HERWIG 


Ps 


t > 8.0 


32.6 


1847 ± 43 


188.5 ± 1.7 


2592 ± 26 


12 


2145 ± 27 


14 


t > 8.5 


24.0 


925 ± 30 


139.0 ± 1.4 


1338 ± 19 


17 


1091 ± 19 


21 


t > 9.0 


16.4 


424 ±21 


95.1 ± 1.2 


630.2 ± 13 


24 


524.1 ± 13 


29 


t > 9.5 


10.1 


179 ± 13 


58.4 ±0.9 


263.8 ±8.3 


36 


229.5 ±8.8 


41 


t > 10.0 


5.5 


76 ± 8.7 


31.8 ±0.7 


105.6 ±5.3 


49 


89.8 ±5.5 


58 


t > 10.5 


2.7 


33 ±5.7 


15.7 ±0.5 


35.1 ±3.0 


73 


35.1 ±3.4 


73 



Table 1: Numbers of events within instanton enhancing cuts chosen such that 
a fraction rj of the QCDINS sample within fiducial cuts (see Sect. 15. g)) is kept. 
Statistical errors are given. The separation power Ps is also shown (see Sect.^. 
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Figure 1: Kinematics of instanton-induced ep collisions. 
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Figure 2: Distribution of the kinematic variable Q'da describing the hard sub- 
process of instanton-induced events. Statistical error bars are given. 
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Figure 3: Distributions of variables calculated from the instanton region, shown 
for the fiducial sample: (a) the transverse momentum, p'^^ , of the current jet, (b) 
EFO multiplicity, N-^-poj (c) EFO track multiplicity, N-^-pT, (d) circularity, C, in 
the hems. Statistical errors are given. 
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Figure 4: Distributions of two of the shape variables for the fiducial sample (a) 
sphericity S, (b) e' . 
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Figure 5: Fisher discriminant t calculated from the variables S, C, p!^^ , N^fOj 
N-EFT and e' . Shown are (a) linear and (b) logarithmic plots for the fiducial sam- 
ple with an additional cut Q'\)p^ < 250 GeV^. In (b), the QCDINS distribution, 
normalised to the predicted fraction of 0.79%, is also shown. 
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Figure 6: Plot of the ratio of the cross sections (Trf/cxth vs. rj, the fraction of 
instanton-induced events remaining after the specific cut in t. The cross sections 
ad and ath refer to the measured cross-section and to the cross-section of instanton- 
induced events according to the QCDINS Monte Carlo generator, respectively. In- 
ner error bars correspond to statistical errors, outer error bars show statistical and 
systematic errors added in quadrature. The scale on the right hand side gives the 
measured cross section in the kinematic region of the fiducial sample. The dashed 
line is a conservative 95% c.l. (see text), corresponding to vj = 4.6%. 
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